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Abstract — It  is  well  known  that,  to  accurately  retrieve  the 
spectrum  of  the  water-leaving  radiance  and  derive  the  ocean 
color  products  from  satellite  sensors,  a  vicarious  calibration 
procedure,  which  performs  sensor  in-orbit  calibration  for  a 
whole  system  (the  sensor  and  algorithms)  is  necessary.  Both 
Sea-viewing  Wide  Field-of-view  Sensor  (SeaWiFS)  and  Mod¬ 
erate  Resolution  Imaging  Spectroradiometer  (MODIS)  have 
employed  in-orbit  vicarious  calibration  procedure  that  uses 
the  in  situ  measurements  with  the  Marine  Optical  Buoy 
(MOBY)  in  the  waters  off  Hawaii.  Such  method  can  also  be 
applied  to  vicarious  inter-calibrate  other  sensors.  In  addition 
to  the  ocean  color  products,  aerosol  optical  property  data  over 
ocean  are  routinely  retrieved  from  both  SeaWiFS  and 
MODIS  measurements.  The  aerosol  retrieval  algorithm  uses 
radiances  measured  at  two  near-infrared  (NIR)  wavelengths, 
at  which  the  ocean  appears  black  due  to  strong  absorption  by 
water,  to  estimate  the  aerosol  optical  properties  and  extrapo¬ 
late  these  into  the  visible.  The  spectral  information  from  two- 
band  measurements  is  used  to  retrieve  the  most  appropriate 
aerosol  models.  With  the  derived  aerosol  models,  the  aerosol 
optical  thickness  can  then  be  estimated  using  the  measured 
signal  at  865  nm.  In  this  paper,  I  outline  the  procedure  for 
the  in-orbit  sensor  vicarious  calibration  for  the  ocean  color 
and  aerosol  products.  Simulations  that  demonstrate  the 
effectiveness  of  the  vicarious  calibration  method  on  the 
derived  ocean  color  and  aerosol  products  are  presented  and 
discussed.  Results  of  sensitivity  studies  that  show  effects  of 
the  calibration  error  at  865  nm,  appropriateness  of  aerosol 
models,  and  the  solar-sensor  viewing  geometry  on  the  accu¬ 
racy  of  the  retrieved  ocean  color  and  aerosol  optical  proper¬ 
ties  are  presented. 

I.  INTRODUCTION 

In  the  ocean  color  remote  sensing,  <10%  of  the  sensor- 
measured  radiance  in  the  visible  at  satellite  altitude  is  from 
the  ocean  near-surface  water,  and  scattering  from  the 
atmosphere  and  the  sea  surface  comprises  the  rest  [1,  2], 
This  places  very  stringent  requirements  on  the  sensor’s 
radiometric  calibration,  particularly  in  the  blue  bands. 
Even  for  a  perfect  atmospheric  correction  (i.e.,  effects  of 
atmosphere  and  surface  are  accurately  removed),  the  rela¬ 
tive  error  in  the  derived  ocean  color  spectrum  will  be  at 
least  10  times  the  relative  error  in  the  sensor  calibration. 
Thus,  for  an  uncertainty  of  5%  in  the  retrieved  water¬ 
leaving  radiance  in  the  blue,  the  uncertainty  in  the  radio- 
metric  calibration  of  the  sensor  must  be  within  -0.5%.  It  is 
impossible  to  achieve  such  an  uncertainty  purely  through 
the  pre-launch  sensor  radiometric  calibration  in  the  labo¬ 
ratory.  Even  if  this  were  possible,  the  launch  process  would 


likely  cause  unknown  changes  in  the  sensor’s  radiometric 
response.  The  key  calibration  procedure  for  the  ocean  color 
mission  is  the  in-orbit  vicarious  calibration  [3,  4].  The 
basic  strategy  is  to  account  for  all  of  the  components  of  the 
top  of  atmosphere  (TOA)  radiance  reflected  from  the 
ocean-atmosphere  system  by  direct  measurement  or  by 
prediction  based  on  in  situ  measurements,  radiative  trans¬ 
fer,  and  scattering  theory.  The  computed  TOA  radiances 
are  then  compared  with  the  results  from  sensor-measured 
radiances.  Any  difference  between  the  measured  and  com¬ 
puted  TOA  radiance  is  attributed  to  error  in  the  calibration 
of  the  sensor,  and  calibration  gain  coefficients  can  be  there¬ 
fore  derived  to  force  the  measurements  and  computations 
into  confluence  (the  computed  TOA  radiance  is  considered 
as  correct  value). 

Gordon  [4]  provided  a  strategy  of  the  in-orbit  vicarious 
calibration  for  the  ocean  color  sensors.  In  this  scheme,  it  is 
assumed  that  the  longest  wavelength  band  in  the  near- 
infrared  (NIR)  is  perfectly  calibrated,  and  then  all  other 
bands  (shorter  than  the  longest  wavelength)  are  in  effect 
calibrated  with  respect  to  this  band.  This  is  the  relative 
calibration.  Wang  and  Gordon  [5]  show  that  as  long  as  the 
calibration  error  at  the  longest  NIR  band  (i.e.,  865  nm)  is 
within  -5%  in  magnitude,  the  in-orbit  vicarious  calibration 
can  produce  the  TOA  radiances  that  are  sufficiently  accu¬ 
rate  to  retrieve  water-leaving  radiances  with  very  good 
accuracy.  This  is  completely  independent  of  the  initial  pre¬ 
launch  calibration  uncertainty  for  all  the  shorter  wave¬ 
lengths  (i.e.,  shorter  than  865  nm).  The  in-orbit  vicarious 
calibration  has  been  successfully  used  to  both  the  Sea¬ 
viewing  Wide  Field-of-view  Sensor  (SeaWiFS)  [6]  and  the 
Moderate  Resolution  Imaging  Spectroradiometer  (MODIS) 
using  the  Marine  Optical  Buoy  (MOBY)  [7]  for  the  ocean 
color  products.  The  in-orbit  vicarious  calibration  scheme 
can  be  also  used  for  the  vicarious  inter-calibration  between 
sensors  [8]. 

In  deriving  the  ocean  color  products  from  ocean  color 
sensors  (e.g.,  SeaWiFS),  the  aerosol  effects  must  be  accu¬ 
rately  estimated  and  removed.  Therefore,  the  SeaWiFS  can 
also  produce  the  aerosol  optical  property  data,  in  particular, 
the  aerosol  optical  thickness  and  Angstrom  exponent,  over 
global  ocean  [1,9]. 

In  this  paper,  I  briefly  describe  the  procedure  for  the  in- 
orbit  sensor  vicarious  calibration  for  the  ocean  color  and 


US  Government  work  not  protected  by  US  Copyright 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

25  JUL  2005 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

In-Orbit  Vicarious  Calibration  for  Ocean  Color  and  Aerosol  Products 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

NOAA  National  Environmental  Satellite,  Data,  and  Information  Service 
Office  of  Research  and  Applications  E/RA3,  Room  102,  5200  Auth  Road, 
Camp  Springs,  MD  20746,  USA 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

1 1 .  SPONSOR/MONITOR' S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM001850,  2005  IEEE  International  Geoscience  and  Remote  Sensing  Symposium  Proceedings 
(25th)  (IGARSS  2005)  Held  in  Seoul,  Korea  on  25-29  July  2005. ,  The  original  document  contains  color 
images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

uu 

18.  NUMBER 
OF  PAGES 

4 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Figure  1.  The  residual  calibration  error  a'(x)  (%)  (panels  (a)  and  (c))  and  ratio  a'(X)/a'(865)  (panels  (b)  and  (d))  as  a  function 
of  the  wavelength  for  a  (865)  values  of  ±5%,  +10%,  and  +15%. 


aerosol  products.  Discussions  that  explain  physics  of  the 
relative  vicarious  calibration  are  provided.  Examples  from 
simulations  that  demonstrate  the  effectiveness  of  the 
vicarious  calibration  method  on  the  derived  ocean  color 
and  aerosol  products  are  presented  and  discussed.  In 
addition,  results  of  sensitivity  studies  that  show  effects  of 
the  calibration  error  at  865  nm  on  the  accuracy  of  the 
retrieved  ocean  color  and  aerosol  optical  properties  are 
presented. 

II.  VICARIOUS  CALIBRATION  PROCEDURE 

By  using  the  reflectance  p=nL/F0  cosOQ  ,  where  L  is  the 
radiance  in  the  given  viewing  direction,  F0  is  the 
extraterrestrial  solar  irradiance,  and  0O  is  the  solar  zenith 
angle,  the  TOA  reflectance  for  the  ocean-atmosphere 
system,  measured  at  a  wavelength  A,  can  be  written  as  [1] 

P,  (A)  =  p,.  (A)  +  pA  (A)  +  t  (A)  pw  (A) ,  (1) 

where  pr  (A) ,  pA  (A) ,  and  pM,( A)  are  the  reflectance 
contributions  from  the  molecules  (Rayleigh  scattering), 
aerosols  (including  Rayleigh-aerosol  interactions),  and 
ocean  waters  (water-leaving  reflectance),  respectively. 
The  t( X)  is  the  atmospheric  diffuse  transmittance  that 
accounts  for  the  effects  of  propagating  water-leaving 
reflectance  from  the  sea  surface  to  the  TOA  [10].  It  is 


computed  with  assumption  that  the  backscattered  radiance 
just  beneath  the  sea  surface  is  uniform,  i.e.,  the  water  body 
is  a  Lambertian  reflector  observed  just  beneath  the  ocean 
surface.  In  the  above  equation,  the  surface  sun  glitter  and 
whitecap  terms  have  been  ignored. 

A  set  of  realistic  aerosol  models  is  needed  for  com¬ 
puting  pA(A)  in  Eq.  (1)  (for  the  atmospheric  correction). 
The  current  SeaWiFS  and  MODIS  ocean  color  data  proc¬ 
essing  uses  12  aerosol  models  for  generating  the  aerosol 
lookup  tables  [9].  They  are  the  Oceanic  model  with  the 
relative  humidity  (RH)  of  99%  (denoted  as  099),  the 
Maritime  model  with  RH  of  50%,  70%,  90%,  and  99% 
(denoted  as  M50,  M70,  M90,  and  M99),  the  Coastal  model 
with  RH  of  50%,  70%,  90%,  and  99%  (denoted  as  C50, 
C70,  C90,  and  C99),  and  the  Tropospheric  model  with  RH 
of  50%,  90%,  and  99%  (denoted  as  T50,  T90,  and  T99). 
These  models  are  from  Shettle  and  Fenn  [11],  and  repre¬ 
sented  well  for  the  maritime  ocean  environment  [12]. 

In  the  vicarious  calibration  (VC),  the  quantities  of  each 
term  on  the  right  hand  side  of  Eq.  (1)  are  estimated  and 
summed,  and  the  sensor  calibration  is  adjusted  to  yield  a 
value  of  pt  (A)  that  is  in  agreement  with  the  estimation. 
Specifically,  the  Rayleigh  reflectance  pr  (A)  is  computed 
using  atmospheric  pressure  and  surface  wind  speed.  The 
water-leaving  reflectance  pw( A)  is  measured  directly  at  the 
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Figure  2.  Errors  in  the  derived  water-leaving  reflectance  and  aerosol  optical  thickness  using  the  Gordon-Wang  algorithm 
with  the  TOA  reflectances  modified  by  the  VC  procedure  for  various  calibration  uncertainties  at  865  nm.  (a)  and  (b)  are 
errors  in  the  derived  water-leaving  reflectance  for  M80  and  T80  models,  and  (c)  and  (d)  are  errors  (%)  in  the  derived  aerosol 
optical  thickness  at  865  nm. 


VC  site  (e.g.,  MOBY)  [7],  pA  (A)  and  t(X)  are  estimated 
using  the  radiative  transfer  theory  with  an  aerosol  model. 
The  VC  procedure  can  be  schematically  outlined  as: 

P,(VC>  (A)  =  ft- (A)  +  pA( A)  +  tpjA) .  (2) 

Computed  VC  aerosol  model  in  situ 

The  above  is  specifically  for  a  VC  site.  We  can  write  the 
sensor-measured  TOA  reflectance  as: 

pJ"1(A)  =  [l+a(A)]pf)(A),  (3) 

where  p\M]  (A)  is  the  reflectance  that  would  be  measured 
by  a  sensor  and  pjn  (A)  is  the  TOA  “true”  reflectance. 
a( A)  is  the  sensor  calibration  errors  (before  VC).  After  the 
VC,  the  sensor-measured  TOA  reflectance  is  modified  as: 

P|C)(A)  =  [l+a'(A)]p*T)(A) ,  (4) 

where  p)c)  (A)  is  the  “Corrected”  sensor-measured  TOA 
reflectance  and  a'{ A)  is  the  residual  calibration  error  after 
VC.  Wang  and  Gordon  [5]  show  that  the  residual 
calibration  error  a'( X)  depends  only  on  the  solar-sensor 
geometry  in  VC  site,  the  VC  aerosol  model,  the  VC  in  situ 
water-leaving  reflectance  measurements,  and  a  (865) 


value.  The  a'(A)  does  not  depend  on  a{X)  for  X  <  865  nm. 
Note  that  a'(865)  =  a(865)  as  it  is  assumed  that  865  nm 
band  is  error  free  in  the  VC  procedure  [4,  5], 

Figure  1  shows  results  of  a'(A)  and  a'(X)/a'( 865)  as  a 
function  of  wavelength  for  a (865)  values  of  ±5%,  +10%, 
and  +15%.  Figs.  1(a)  and  1(b)  are  results  of  a'(X)  and 
fl'(A)/a'(865)  from  a  VC  using  a  correct  aerosol  model 
(model  099),  while  Figs.  1(c)  and  1(d)  are  results  from  an 
incorrect  VC  aerosol  model  (true  099,  but  M99  for  VC). 
In  Fig.  1,  the  VC  has  been  carried  out  specifically  for 
aerosol  optical  thickness  of  0.05  with  “true”  aerosol  model 
of  099,  and  solar-sensor  geometry  of  20°,  20°,  and  90°  for 
the  solar-zenith,  sensor-zenith,  and  relative  azimuth  angle, 
respectively.  Figs.  1(b)  and  1(c)  also  plot  the  curves  of 
(A/865)  and  (a/865)  that  demonstrate  basic  physics  of 
the  VC  procedure — the  ratio  of  residual  calibration  error 
fl'(A)/a'(865)  varies  according  to  more  or  less  the  inverse 
of  the  Rayleigh  scattering,  i.e.,  (A/865 )  .  The  a'(A)  in  the 
visible  is  about  an  order  less  than  a'( A)  at  865  nm.  Thus, 
using  the  outlined  in-orbit  VC  scheme,  for  calibration 
uncertainty  of  ~5%  at  the  865  nm,  the  calibration 
uncertainty  in  the  blue  can  be  reduced  to  within  -0.5%. 


III.  EFFECTS  ON  THE  DERIVED  OCEAN  COLOR 
AND  AEROSOL  PRODUCTS 

To  assess  the  efficacy  of  the  VC  scheme,  simulations 
have  been  carried  to  evaluate  the  accuracy  of  the  derived 
water-leaving  reflectance  pw(A)  and  aerosol  optical 
thickness  to  (A)  using  the  TO  A  reflectances  that  have  been 
modified  with  the  VC  gains,  i.e.,  p)c)  (A)  in  Eq.  (4),  for  the 
various  VC  cases  (e.g.,  VC  aerosol  models,  solar-sensor 
geometry,  etc.),  and  various  retrieval  cases  (e.g.,  aerosol 
optical  properties,  solar-sensor  geometries,  etc.).  For  the 
VC  procedure,  it  was  assumed  that  an  incorrect  aerosol 
model  M99  was  used  for  a  “true”  099  model  (realistic 
case).  This  is  the  same  case  as  in  Figs.  1(c)  and  1(d).  With 
the  derived  a'{  A)  forgiven  a(865)  values  (through  the  VC 
procedure),  the  TOA  reflectance  p'C)  (A)  can  be  computed 
according  to  Eq.  (4).  In  generating  the  testing  data 
(A) ,  the  Maritime  and  Tropospheric  models  with  RH 
of  80%  (M80  and  T80)  [11]  were  used.  These  two  models 
are  different  from  the  12  aerosol  models  used  in  generating 
aerosol  lookup  tables  for  the  atmospheric  correction  and 
aerosol  retrievals.  With  the  generated  p)c)  (A)  values,  the 
Gordon  and  Wang  [1]  atmospheric  correction  and  aerosol 
retrieval  algorithm  was  performed  to  derive  the  water¬ 
leaving  reflectance  pw( A)  and  aerosol  optical  thickness  at 
865  nm  xa  (865) . 

Figure  2  provides  examples  of  the  retrieval  results. 
Figs.  2(a)  and  2(b)  are  uncertainties  in  the  derived  water¬ 
leaving  reflectance  at  443  nm  for  the  M80  and  T80  models, 
while  Figs.  2(c)  and  2(d)  are  errors  (%)  in  the  retrieved 
aerosol  optical  thickness  for  M80  and  T80  models.  These 
are  for  case  of  aerosol  optical  thickness  ra  (865)  of  0.1  and 
for  a  specific  sensor-viewing  geometry  (sensor-zenith  angle 
of  46°  and  relative  azimuth  angle  of  90°)  with  solar  angles 
varying  from  0°-80°.  They  are  corresponding  to  the  VC 
cases  for  a  (865)  values  of  0%,  ±5%,  ±10%,  and  +15%, 
respectively. 

Results  in  Figs.  2(a)  and  2(b)  show  that  the  water¬ 
leaving  reflectance  at  443  nm  can  be  retrieved  accurately  as 
long  as  the  sensor  calibration  uncertainty  at  865  nm  a  (865) 
is  within  ~5%.  In  fact,  for  the  M80  model,  errors  in  the 
derived  water-leaving  reflectance  are  all  within  0.002  and 
mostly  within  ±0.001  for  a(865)  up  to  ±15%.  For  the  T80 
model,  except  some  cases  with  larger  solar  zenith  angles 
for  a(865)  =  +15%,  the  derived  water-leaving  reflectances 
are  all  accurate  within  ±0.002. 

On  the  other  hand.  Figs.  2(c)  and  2(d)  show  that  error  in 
the  derived  aerosol  optical  thickness  at  865  nm  is  almost 
twice  the  magnitude  as  the  calibration  uncertainty  a  (865), 
e.g.,  5%  in  a(865)  leading  to  ~10%  error  in  the  derived 
aerosol  optical  thickness  xa  (865) .  This  is  because  the 
aerosol  scattering  contributes  -50%  of  the  sensor-measured 
TOA  reflectance  (Rayleigh  scattering  contributes  another 
50%).  Calibration  error  in  the  TOA  reflectance  at  865  nm 
is  equivalent  to  double  the  effect  on  the  aerosol  component. 


Therefore,  to  retrieve  the  aerosol  optical  thickness  within 

error  of  10%,  a (865)  has  to  be  within  5%. 
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